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SUMMARY 

A simple engineering correlation scheme Is developed to predict the vari- 
able property effects on dilute species laminar forced convection mass transfer 
applicable to all vapor molecules or Brownian diffusing small particles, 
covering the surface to mainstream temperature ratio of 0.25 < T w /T e 1 4. 

The accuracy of the correlation Is checked against rigorous numerical forced 
convection laminar boundary layer calculations of flat plate and stagnation 
point flows of air containing trace species of Na, NaCl, NaOH, Na 2 S 04 , K, KC1, 
KOH, or K 2 SO 4 vapor species or their clusters. For the cases reported here 
the correlation had an average absolute error of only 2 percent (maximum 13 
percent) as compared to an average absolute error of 18 percent (maximum 54 
percent) one would have made by using the "constant-property" results. 


INTRODUCTION 

Although considerable temperature variations are encountered across the 
boundary layer (BL) In many heat and mass transfer situations, no simple pro- 
cedure to account for variable (nonconstant) property (ncp) effects on mass 
transfer rate exists In the literature. For such cases, the application of 
the constant property (cp) analytic solutions, or the experimental data, may 
lead to unacceptable errors. Readily available corrections are urgently needed 
to Improve the accuracy of transfer rate predictions and to facilitate engi- 
neering design. 
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Most of the nonconstant property results for heat transfer prediction 
purposes Indicate that very simple corrections generally suffice over a moder- 
ate temperature range. A thorough review of the current literature as well as 
a detailed discussion of the subject Is given In reference 1 both for liquids 
and gases for laminar and turbulent flows. One of the common schemes used to 
correct for variable properties Is the property ratio scheme. For gases, the 
temperature-dependent-property effects for heat transfer can usually be ade- 
quately correlated by: 


_ ( ‘ q w ) 
ncp.h = 

v M w'cp 


Nu. 


Nu 


h,cp 



h,cp 


( 1 ) 


where the constant-property quantities are evaluated at free-stream temperature 
for external flows. 


For mass transfer across boundary layers with steep temperature gradients, 
however, a similar correlation scheme Is lacking, and the conventional heat and 
mass transfer analogy Is usually and unjustifiably employed. Many engineering 
fields Involving mass transfer such as CVD, surface catalyzed combustion, salt 
deposition on gas turbine blades associated with hot corrosion, etc. deal with 
nonlsothermal BLs. A carefully formulated mass transfer theory, therefore, 
explicitly necessitates tractable variable-property correlations to facilitate 
efficient and accurate mass transfer rate predictions (2 and 3). 

A correlation of variable property effects on mass transfer for 
chemically-frozen laminar BLs has been proposed In reference 4 using an equiv- 
alent "blowing (suction) parameter" for a number of vapor species. The present 
paper takes a much simpler approach similar to the property ratio scheme as 
exampllf led by equation 1 and extends the applicability to larger particle 
sizes Including vapor species (10-3 < Le < 1). 


NOMENCLATURE 

a power of temperature ratio for variable property effects on heat transfer 
equation (1) 

b power of temperature ratio for variable property effects on mass transfer, 
equation (14) 

c p specific heat at constant pressure per unit mass of mixture 
0 Brownian (Flck) diffusion coefficient 
Eu Euler number, d In u e /d In x 

F function describing variable property effects, equations (1) and (7) 

j" mass (diffusion) flux 

Kn Knudsen number based on particle diameter 
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k Boltzmann constant 

Le Lewis number (ratio of Brownian (Flck) dlffuslvlty to carrier gas thermal 
dlffuslvlty) 

Nu Nusselt number 

n power of temperature for Brownian (Flck) diffusion coefficient 

Pr Prandtl number (ratio of carrier gas kinematic viscosity to thermal 
dlffuslvlty) 

q" heat flux 

Re Reynolds number based on x 

r power of temperature for thermal conductivity 

Sc Schmidt number (ratio of gas kinematic viscosity to particle Brownian 

(Flck) dlffuslvlty) 

St Stanton number 

T absolute temperature 

u velocity along x 

W dimensionless mass fraction of particles In prevailing mixture, 
equation (3) 

x distance along the surface 

y distance normal to surface 

a power of temperature for specific heat 

(1 power of temperature for viscosity 

g Lennard- Jones molecular Interaction energy (well-depth) parameter 
n dimensionless BL coordinate, equation (4) 

e dimensionless temperature, equation (2) 

X thermal conductivity of carrier gas mixture 

y viscosity of carrier gas mixture 

p density of carrier gas mixture 

w mass fraction of particles In the carrier gas mixture 
Q collision Integral 
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Subscripts: 

av averaged quantity 

cp constant property 

0 pertaining to Brownian (Flck) diffusion coefficient 

e outer edge of BL, mainstream conditions 

h heat transfer 

1 species 1 

j species j 

m mass transfer or outer edge of the mass transfer BL (e.g. T m ) 
ncp nonconstant (variable) properties 
w wall (surface) conditions 

Superscripts : 

' differentiation with respect to n 
* pertaining to dimensionless temperature, kT/e av 

Miscellaneous: 

BL boundary layer 

CVO chemical vapor deposition 


CORRELATION APPROACH 

Here we concentrate on chemically-frozen, gaseous, laminar BLs where the 
host gas thermodynamic properties and the trace gaseous species transport prop- 
erties are assumed to have simple power-law temperature dependencies given by: 


P ~ T 1 , X ~ T r , y ~ T®, Cp ~ T a , D ~ T n 


If one applies the conventional self-similarity transformation to the two- 
dimensional laminar BL equations and uses the nondlmenslonal variables defined 
by: 
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where n Is 


n 



Re 


1/2 
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(4) 


equation (1) can then be reexpressed as: 


ncp,h 



e' 

w.cp 


(5) 


Tables of eigenvalues, e w* are presented for a family of wedge solutions by 
Brown and Oonoughe (ref. 5) when air Is used as the host gas (r = 0.85, 8 = 
0.70, a = 0.19, Pr = 0.70). Based on these solutions, equation (5) combined 
with equation (1) reduces to: 


ncp,h 



( 6 ) 


where the best temperature ratio exponent, a, Is given In table 1 as obtained 
from reference 1 . 


Similar to equation (1), the effect of variable properties on mass trans- 
fer, In the absence of thermal diffusion ( "thermophoresis" for small parti- 
cles), can be described as: 
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Combining with equations (3) and (4) we have: 



which for air reduces to: 


( 8 ) 
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The eigenvalues for mass transfer rate Including variable properties, W^, 
are obtained from the numerical solutions of the BL equations for any size of 
diffusing trace species, as outlined In reference 6. The eigenvalues for heat 
transfer rate, e^, generated as by-products by the calculations of 
reference 6, agree excellently with the values reported In reference 5, and 
confirm the temperature ratio exponent given In table 1. 

Based on equation 8, attention Is directed toward correlating the mass 
transfer eigenvalue ratio, Wj,/Wj, cp . It Is known that In the limiting case 
of Sc » 1 for small particles, tne mass transfer BL Is very thin compared to 
the momentum (or heat) transfer BL. Therefore, constant-property results 
obtained by evaluating the properties at the wall temperature will give correct 
answers. It Is also known from asymptotic laminar BL theory for Sc » 1 that. 
In the absence of thermophoresis, Wj, Is proportional to Scl/3 (ref. 7). 

Hence, for Sc » 1 : 


w; “w.cp,T.T u 

w - w 1 

w,cp w,cp,T=T e 

which for air reduces to: 



Equation (11) Is confirmed also by the calculations of reference 6. 

In the limiting case where Le -» 1 , however, It Is expected that F ncp m 
comes close to F ncp h. the difference being due primarily to the temperature 
dependence of the diffusion coefficient of the species (D ~ T n ) . Note that 
In the case of air as the host gas, the temperature dependence of the Prandtl 
number Is: 


Sc 


Sc 


1/3 

'w 

1/3 



( 10 ) 


p r ~ xB+a-r = j0«04 


( 12 ) 


and the temperature dependence of the Schmidt number Is: 


Sc ~ T&+l-n = t 1 «70-n 


(13) 
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Hence, for Le = 1 , F nC n m = F ncp> h only If n = 1.66. This should be 
compared with n = 1.652'used for’the vapor species considered In reference 4. 
In fact, n = 1.66 Is a good approximation for species diffusing In a gas with 
T* = kT/e av greater than about 4, where e av = (c^gj ) 1 / 2 (see also the 
"OISCUSSION" section). Based on the arguments given above, we propose the 
following simple correlation: 



which for air reduces to: 


W 

w w 


W 

w,cp 



1.70-n 



b-n 


(15) 


where b Is a constant and T m Is the temperature at the outer edge of the 
mass transfer BL. In the absence of thermophoresis T m Is well approximated 
by: 

1/3 

Tm s ^w + *-e av (T e - T w ) (16) 

where Le av Is evaluated at T av = (T w + T m )/2. In equations (14) and 
(15) the first term In the parentheses accounts for the temperature dependence 
of Sc, and the second term In the parentheses accounts for the property vari- 
ation across the mass transfer BL. Note that for Le = 1 , equation ( 8 ) 
becomes : 



which Is Independent of n, and, for air, reduces to: 


/TjA b '’- 85 
r ncp,m (t e ) 


If one assumes that F nC p >m = F nC p ^ for Le = 1 , then the value of b 
Is fixed and Is given by D = a + 1.85, as tabulated In table 2. The accuracy 
of equation (15) with the values of b given In table 2 Is checked against 
rigorous numerical calculations as shown In figure 1 , where air Is the carrier 
gas. Indeed, our suggested correlation Is found to be generally applicable 
and successful for many different vapor species like Na, NaCl, NaOH, Na 2 S 04 , 

K, KC1, KOH, K 2 SO 4 and clusters made of these molecules. The cases considered 
covered the ranges 0.25 < T w /T e < 4 and 10 ~ 3 < Le < 1 for both the flat plate 
(Eu = 0) and two-dimensional stagnation point (Eu = 1) flows. Our correlation 
gave a mean absolute error of only 2 percent (maximum 13 percent) as opposed 
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to a mean absolute error of 18 percent (maximum 54 percent) If the "constant- 
property" results were used. 


DISCUSSION 

The generality of the correlation suggested above Is accomplished by 
condensing the temperature dependence of the Individual transport characteris- 
tics of different species Into a single parameter n. Therefore, In order to 
utilize the correlation one needs to be able to readily obtain the temperature 
dependence of the diffusion coefficient for any species In question. For vapor 
molecules the temperature dependence of the diffusion coefficient Is given by 
the Chapman-Enskog theory as: 


D 


t 3/2 


= T' 


(19) 


where Qq Is the collision Integral which Is a function of the dimensionless 
temperature T* = kTav/eav evaluated at the average mass transfer BL 

temperature T flv = j (T m + T w ). From equation (19) one obtains that: 


n 


3 

2 


d 1n Q p 

d In T* 


( 20 ) 


The calculation of the collision Integral Is cumbersome and Is given In 
reference 8. In order to facilitate the calculation of n, the tabulated data 
given In reference 7 Is numerically differentiated to get an approximation to 
d In Qg/d In T*. A curve, based on the n expression given In reference 8, Is 
fitted through these points as shown In figure 2. The analytical expression 
describing this curve-fit Is given by: 

d In s 

" d In T* s °' 150 + 0,385 exp [~ O* 637 log-, gT* + 0.340)] (21) 

As can be seen from figure 2, for high temperature situations with air as the 
carrier gas where T* Is usually greater than about 4, n = 1.66 Is a good 
approximation for many vapor species. This supports the assumption made ear- 
lier that F nC p >m = Fpcp^h for Le = 1 . 

In the extreme limit where the diffusing particles are large enough to be 
treated as If they were In a continuum (Sc » 1, Kn « 1), the temperature 
dependence of the diffusion coefficient Is given by the Stokes-Elnsten theory 
as: 



( 22 ) 
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and therefore, n = 1 - 13. In this case we have from equations ( 8 ) and (10): 



13+1-1+13 3+513 



which for air reduces to: 


F 


ncp.m 



-1 .08 


(23) 


(24) 


For Intermediate size particles which are not In the continuum regime 
there Is no straightforward method to obtain n. For such cases, the diffusion 
coefficient can be calculated as outlined In reference 9, and then n can be 
obtained from: 


n 


in r°< T J /D <vi 

ln(T„/T. J 


(25) 


One should bear In mind, however, that whenever Sc Is large enough to vali- 
date the Sc » 1 asymptotic theory obtained for constant properties, the mass 
transfer BL becomes extremely thin compared to the momentum (or heat) transfer 
BL. For such cases we suggest that Wy be exactly obtained from the 
asymptotic theory only with properties evaluated at the "wall" conditions. 
Instead of using the variable property correction recommended above. 

The effect of variable properties on mass transfer naturally depends on 
(1) the magnitude of the difference between the mainstream temperature, T e , 
and some average temperature that the mass transfer BL experiences, and (2) 
the sensitivity of the species diffusion coefficient to temperature (namely, 
n). Therefore, similar to the argument given for Sc » 1 (Le « 1) cases, 
for Sc « 1 (Le » 1) most of the mass transfer BL experiences a constant 
temperature profile equal to the mainstream temperature, T e . Considering 
the fact that n Is almost the same for most vapor species, we expect, there- 
fore, that the effect of variable properties on vapor species with Le > 1 
will be even smaller than the effect observed for vapor species with Le < 1 . 
Table 3 lists the exact numerical results obtained for two vapor species, 
namely K 2 SO 4 and Na. At 1500 K K 2 SO 4 has Le = 0.344 (smallest among 
the vapor species tried) and Na has Le = 0.807 (largest among the vapor 
species tried). As expected, the effect of variable properties for the Na 
case with the larger Le Is less than for the K 2 SO 4 case. For the other 
vapor species tried, with Intermediate Le, It Is confirmed that the effect of 
variable properties ranks according to the magnitude of Le. Based on the 
above argument and given the fact that variable property effects for Na are 
already negligibly small, we propose that F nC p m for vapor species with 
Le > 1 should be taken as unity. 

The following Is an outline of the recommended procedure that one should 
follow to obtain the correction factor, F nC p m , for any given nonlsothermal 
laminar forced convection trace species mass ’transfer application. It Is 
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assumed that conventional BL procedures have already been fully exploited to 
determine the mass transfer rate for constant properties, l.e., -jJ CD or 
N u m,cp or St m>C p are accurately known. When air Is the carrier (host) gas: 

Step 1: Calculate the Le of the trace species. If Le > 1, then take 

^ncp.m = 1 • It Le < 10 -3 , then use the constant-property 
solutions with properties evaluated at T w . 

Step 2: If the trace species Is a vapor molecule, determine 

d In Qq /d In T* from equation (21) (or fig. 2) and calculate n 
from equation (20). 

Step 3; If the trace species Is a small particle (cluster), determine n 
from equation (25) . 

Step 4: Calculate T m from equation (16). 

S tep 5: Obtain the proper value of b from table 2. If It Is not a flat 

plate or a stagnation point flow (e.g., wedge flow, etc.) use some 
average value of b from table 2. 

Step 6: Use equation (15) to obtain Wy/wJ, fCp . 

Step 7: Use equation (9) to obtain F ncpm . 

££§P 8; Use equation (7) to obtain the mass flux Jy, Nu m or St m 

corrected for variable property effects. 

If the carrier gas Is not air, the formulation In the paper Is kept gen- eral 
to allow Steps 1 to 8 (except Step 5) to be repeated with the proper values of a, 

0 and r corresponding to the temperature dependencies of carrier gas properties. 
However, the values of b as listed In table 2 (Step 5) are only for air and have 
to be determined for other gases. 


CONCLUDING REMARKS 

Interfaclal mass transport rate formulations Including thermal (Soret) 
diffusion ("thermophoresis" for small particles) effects require explicit 
knowledge of both F ncp k and F ncp>m as shown In references 2 and 3. 

Toward this end, a simple correlation scheme Is developed to predict the vari- 
able property effects on mass transfer. The applicability range of the cor- 
relation covers 0.25 < T w /T e < 4 and 10 -3 < Le < 1 . When checked against 
rigorous numerical calculations for (1) dilute vapor species Na, NaCl, NaOH, 
N f2 30 4* K, KC1 , KOH, K 2 S0 4 , (2) their clusters, and for (1) laminar flat 
plate flows and (2) laminar stagnation point flows, It Is found that the 
correlation had an average absolute error of only 2 percent (maximum 13 per- 
cent) as opposed to an average absolute error of 18 percent (maximum 54 per- 
cent) one would have made by using "constant-property" results. We propose 
that constant property solutions evaluated at the "wall" temperature be used 
for larger particles than those with Le = 10" 3 and expect the variable 
property effect to be negligibly small for cases where Le > 1 . 
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TABLE I. - EXPONENT OF TEMPERATURE RATIO FOR VARIABLE 
PROPERTY EFFECTS ON HEAT TRANSFER, a 




era 

T e >T « 

Flat plate 

Two-dimensional stagnation point 

(Eu = 0) 
(Eu = 1) 

-0.01 

.10 

0 

.07 


TABLE II. - EXPONENT OF TEMPERATURE RATIO FOR VARIABLE 
PROPERTY EFFECTS ON MASS TRANSFER, b 




na 

BB 

Flat plate 

Two-dimensional stagnation point 

(Eu = 0) 
(Eu = 1) 

ID 

m 


TABLE III. - EXACT CALCULATIONS OF VARIABLE PROPERTY EFFECTS 
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Figure 1. - Accuracy of correlations for flat plate and two-dimensional stagnation point flows in the range 
a 25< T w /T e « 4 and 10‘3 £ Le £ 1. 
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Figure 2. -d In Q D /d In T* as a function of T* obtained by numerical differentiation of tabulated Q n 
values of ref. 7. A curve fit through the points is given by eq. (211. u 
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